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Abstract The system of one-dimensional symmetric simple random walks, in which none of walkers have 
met others in a given time period, is called the vicious walker model. It was introduced by Michael Fisher 
and applications of the model to various wetting and melting phenomena were described in his Boltzmann 
medal lecture. In the present report, we explain interesting connections among representation theory, proba- 
bility theory, and random matrix theory using this simple diffusion particle system. Each vicious walk of N 
Qh walkers is represented by an TV-tuple of nonintersecting lattice paths on the spatio-temporal plane. There 

Ph ' is established a simple bijection between nonintersecting lattice paths and semistandard Young tableaux. 

Based on this bijection and some knowledge of symmetric polynomials called the Schur functions, we can give 
a determinantal expression to the partition function of vicious walks, which is regarded as a special case of the 
Karlin-McGregor formula in the probability theory (or the Lindstrom-Gessel-Viennot formula in the 
cnumerative combinatorics). Due to a basic property of Schur function, we can take the diffusion scaling limit 
of the vicious walks and define a noncolliding system of Brownian particles. This diffusion process solves 
the stochastic differential equations with the drift terms acting as the repulsive two-body forces proportional to 
the inverse of distances between particles, and thus it is identified with Dyson's Brownian motion model. 
In other words, the obtained noncolliding system of Brownian particles is equivalent in distribution with the 
, eigenvalue process of a Hermitian matrix-valued process. 

i/-) ! 1 Vicious Walks, Young Tableaux and Schur Functions 

O , 

C"| ■ Let ^{S(t)}t = o.i,2...., P x ) be the TV-dimensional Markov chain starting from 

x = (x±, X2, • ■ ■ , Xn), such that the coordinates Si(t), i — 1, 2, . . . ,N, are independent simple random walks on 
Z. We always take the starting point x from the set 



oo 



a 



Z N 



x 



< 



jx = (xi,X2, • • • ,Xn) E (2Z) N : x\ < x 2 < ■ ■■ < Xn\- 



We consider the condition that any of walkers does not meet other walkers up to time T > 0, i.e. 

S^t) < S 2 {t) < ■ ■ ■ <S N (t), t = l,2,...,T. (1.1) 

We denote by Q T the conditional probability of P x under the event At = |<Si(£) < <?2(i) < • • • < Sjv(t), t = 

o,i,. ..,r}. M. Fisher called the process ^{S(£)} t =o,i,2,...,T, the vicious walker model in his Boltzmann 

medal lecture [4]. 

We will assume the initial positions as 

$(0)=2(i-l), i = l,2,---,N (1.2) 
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Figure 1: An example of vicious walk in the case Y = 4, T = 6. 



in Sections 1 in this report. 

Each realization of vicious walk is represented by an TV-tuple of nonintersecting lattice paths on the 
1+1 spatio-temporal plane, Z x {0, 1, ■ ■ ■ , T}. An example is given by Figure 1 in the case that four walkers 
(N = 4) perform a noncolliding walk up to time T = 6. 
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Figure 2: (a) Young diagram and (b) Young tableau T 
corresponding to the vicious walk in Figure 1 . 



Figure 3: A pair of conjugate YD's. 



Bijection between such nonintersecting lattice paths and semistandard Young tableaux (SSYT), T = 

(T(i,j)), is established by the following procedure [8, 17]. 

(1) For 1 < j < N, let 

Lj = the number of leftward steps among T steps of the j-th walker. 

Draw a collection of boxes with N columns, in which the number of boxes in the j-th column is Lj. (We 
number columns from the left to the right.) Since L = (L\, L2, £3, L4) = (3, 3, 2, 1) in the walk shown in 
Figure 1 , we draw the collection of boxes as shown in Figure 2 (a) for this example. 

(2) For each walker, we label each leftward step by the integer £ {1, 2, • • • , T}, which is the time when that 
leftward step was done. See Figure 1, in which labels of leftward steps are indicated by integers in small 
circles associated with the line segments showing leftward steps. Then for the j-th column of the collection 
of boxes, fill the boxes by the labels of leftward steps of the j-th walker, from the top to the bottom, 
1 < j < N. For the walk given in Figure 1, we have the boxes with integers shown in Figure 2 (b). Let 

T(i,j) — the integer in the box located in the z-th row and j-th column. 

For example, T(l, 3) = 4 and 7~(3, 1) = 5 in this case. 
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Remark 1. The above procedure with the nonintersecting condition (1.1) guarantees the inequalities 

L 1 >L 2 >--->L N , (1.3) 

and 

T(i,j) < T(i + strictly increasing in each column, 

^"(*>j) < + 1), weakly increasing in each row. (1.4) 

Assume that the number of rows in the collection of boxes is £. Let 

Aj = the number of boxes in the i-th row, i = 1, 2, • • • , £ . 

Then the inequalities (1.3) imply 

Ai > A 2 > ••• > Xe- (1.5) 

The collections of boxes with such conditions concerning the numbers of boxes in rows (1.5) (and in 
columns (1.3)) are called Young diagram (YD). The number of rows I in the YD is called the length 
of YD. In the present situation, £ < T, in general. (In our example in Figure 2 (b), I = 3 for T = 6.) 
YD's with integers with the conditions (1.4) are called semistandard Young tableaux (SSYT). 

Remark 2. YD with A^ boxes in the i-th row, 1 < i < £, is said to be the YD of shape A = (Ai, • • • , A^). 
The YD with the shape L = (Li, ■ • • , L N ) is regarded as the conjugate of the YD with the shape A and 
denoted by 

L = A. 

As shown in Figure 3, they are mirror images with respect to the diagonal line. 
A sequence of integers with the condition (1.5), that is, 

A = (Ai, A 2 , • • • , A T ) with Ai e N = ja; e Z, x > oj, 1 < i < T, Ai > A 2 > • • • > A T , 

is regarded as a partition of an integer n = J2i=i We introduce a set of T variables z = {z\, z 2 , • • • , Zt} € C T 
and define a monomial 

z r = n zt (i.j) 

T 

n# of times that the integer k occurs in T 
z k 

k=l 

For example, the monomial corresponding to the SSYT T shown in Figure 2 (b) is 

z r = z 2 x z 3 x z 4 x z 6 

X Z4 X Z4 X Zq 

x z 5 x zq 

3 3 
= Z2Z3Z 4 Z§Z 6 . 

Notice that for one YD with a given shape A, there are different ways of filling boxes with integers to make 
SSYT's satisfying the conditions (1.4). For each YD with shape A, we define a polynomial of z = {zi, • • • , zt} 
by summing z r over all SSYT defined on the YD: 

s\{zi,z 2 ,---,z T ) = zT - 

T:all SSYT with the same shape A 

This polynomial is called the Schur function indexed by (the partition/ YD with shape) A on z = (z\, ■ ■ ■ , zt)- 
We can prove the two formulae (Jacobi-Trudi formulae). The first one is the following. 
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Lemma 1 



sa(zi, ■ • • ,Zt) = 



det 

l<i,3<T 



Xj+T-j 



det 

l<i,j<T 



where the denominator is the Vandermonde determinant evaluated as the product of differences, 



det 

l<i,3<T 



T-j 



n ~ z ^- 

l<i<j<T 



(1.6) 



This formula clarifies that the Schur functions are symmetric polynomials in z = (z\, ■ ■ ■ , Zt)- 

For the second formula, we define the polynomials ej(zi, ■ ■ ■ , Zt)'s as the coefficients in the expansion 

T T 
3=0 



(1.7) 



Then 



ej(zi,---,zr) 



i=i 



sum of all monomials in the form Zi 2 ■ ■ ■ zu 



for all strictly increasing sequences 1 < i\ < i 2 < ■ ■ ■ < ij < T. 



ej(zi, ■ ■ ■ , Zt)'s are also symmetric polynomials in z\,---,Zt and called the j-th elementary symmetric 
polynomials. 















Lemma 2 Assume that the conjugate 


ofX 


is given by 


A = (Ai 


■ ■ ■ , \m) with length N. Then 


s\(zi, ■ ■ ■ 


,zt) 


= det 

l<i,j<N 


e \j+(i-j 


} (zi, ■ ■ ■ , z T ) 

















More details for YD, SSYT and symmetric polynomials, see e.g. Fulton (1997) [5]. 

Now we go back to the vicious walker model. We notice a simple relation between the partition L = 
(Li, ■ ■ ■ ,Ln) and the final positions of the N vicious walkers at time T, y = (yi, • • • , Un) = (Si(T), • • • , Sn(T)), 
given by yi = T — 2Li + 2{i — 1), 1 < i < N, on the initial condition (1.2). Then we have established the 
following relation between the vicious walks and YD/SSYT/Schur functions. 



Set the initial positions as x = (0, 2, • • • , 2(N - 1)). For y such that y e if T e 2N, y + 1 = 

(2/i +1, ■ ■ ■ ,y N + 1) e Z^, if T + 1 e 2N, let L = (Li, ■ • • , L N ) with L 4 = (T + x 0i - y 4 )/2, 1 < i < N, and 
A = L. Then 

• positions y of vicious walkers 

at the finial time T <= 

• a realization of vicious walk 
from S(0) = x to S(T) = y <*= 

• a set of all vicious walks 

from S(0) = x to S(T) = y <*= 



YD with the shape A 

an SSYT T with the shape A 

a Schur function s\(zi, • • • , zt). 



For y = (yi, y 2 , • • • , y N ) such that y e Z^, if T e 2N, y+1 e Z^, if T+l e 2N, and x = (0, 2, • • • , 2(iV-l)), 
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define the number 

Mn,t{y) = total number of distinct realizations of vicious walk of N walkers 
from the positions S(0) = xo to S(T) = y. 
The above relations prove the following identity. 



Assume that L, = TjL^L— ^ — — ,L = (L\, ■ ■ ■ , £jv), and A = L. Then 


Mjv,T(y) = s\(z lr --,z T ) 






Zl=Z 2 = --- = ZT = l 



Lemma 1 gives the following estimate for Mjv^y), via appropriate q- factorization and the formula for 
the Vandermonde determinant (1.6). 



Proposition 3 








M N . T {y) 


= lim s x {l,q,q 2 , ■ ■ 

qr->l 


,9 T - 1 ) 






= limoELiCfc- 1 )^ 

9-1 


n - 

l<i<j<T 


.-Aj+i-i _ ! _ Xi-Xj+j-i 
qj-i-l 11 j-i 

H l<Kj<T J 



On the other hand, by definition (1.7), it is easy to see that 



ej(z!,---,z T ) 
Then Lemma 2 gives the following. 



z\ — — zt — z 



jKT-jy. 



z T . 



Proposition 4 



Mjv,T(y) = 



det 

l<i,j<N 

det 

l<i,j<N 









y 


= det 




l<i,j<N 





Xj+i- jj 



{T + 2(i-l)- yj }/2 



2 Determinantal Formula for Nonintersecting Paths 

Since we have assumed the initial positions as (1.2), we can see that the (i, j)-element of the matrix in the 
determinant in Proposition 4 is 

rm , | = #1 lattice path from 2(i — 1) at time to w, at time t\ 

{T + 2(i-l)- yj }/2j ^\ K ' y] J 

= w(path) 

all lattice paths: (2(i - 1), 0) ~> (j/j , T) tu(path) = l 

If we define an appropriate weight function w(path) on single lattice paths, the summation of w(path) will give 
the Green function of single lattice paths, 

G((x,0),(y,T)) = Yl tu(path). 

all lattice paths: (x,0) ~> (y,T) 
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Proposition 4 can be regarded as a special case of the Karlin-McGregor formula in the probability theory 
[10, 11], and the Lindstrom-Gessel-Viennot formula in the enumerative combinatorics (see [19, 17] and 
references therein). In order to explain this fact, here we introduce some definitions and notations for describing 
lattice paths. 

Let V = {vertex}, E = {directed edge}, D = (V, E) = an acyclic directed graph, where acyclic means that 
any cycles of directed edges are forbidden. For u, v e V, 

a lattice path u — > v = a sequence of directed edges from u to v, 
V(u, v) = the set of all lattice paths from u to v. 

A weight function w : E — > Z[[x e : e € E]] is introduced, where Z[[x e : e € E\] denotes a ring of formal power 
series of {x e : e e E} and the weight on a lattice path P is defined by w(P) = J^J w(e). Then the Green 

eeP 

function of lattice paths from u to v is defined by 

G(u,v)= Yl 

Let / = {ui, U2, ■ ■ ■ , uat}, J = {vi, V2, ■ ■ ■ , vn} with Ui,Vi € V, i = 1,2, •••,7V. The sets I, J are ordered; 
mi < U2 < ■ ■ ■ < un, vi < V2 < ■ ■ ■ < vn- Then we consider a set of iV-tuples of lattice paths 

V(I, J) = jp = (Pi, ■ • • , P N ) : Pi e 7>(«i, «0, i = 1, 2, • • • , N}. 

The weight for each iV-tuple of lattice paths is given by 

N N 

»(p)=n^)=nn ^ 

i=l i=l eGPi 

We say 'lattice paths P and Q intersect', if P and Q share at least one common vertex. Then, for ordered 
sets of vertices I and J, we say 'I is D-compatible with J' in the case that, whenever u, < itj in / and 
«i < t)j in J, every lattice path P € V(ui,vj) intersects every lattice path Q e V(uj,Vi). A set of TV-tuples of 
nonintersecting lattice paths is denoted by 

Vo(I, J) = |p € P(/, J) : any lattice paths in P do not intersect with others j, 

and the Green function of TV-tuples of nonintersecting lattice paths is defined by 

G mmiBt {i,j)= ™( p )- 

Pev (i,J) 



Theorem 5 Let I = (ui, ■ ■ ■ , it/v) and J — {vi, • • • ,ujv) be two ordered sets of vertices in an acyclic graph 
D. If I is D-compatible with J, then the Green function of the nonintersecting N -tuples of lattice paths is 
given by 



Gnonint(P J) = det 

l<i,j<N 



G(ui,Vj) 

where G(u, v) denotes the Green function of single lattice paths from u to v on D. 



The proof of Theorem 5 is given in Appendix A following Stcmbridgc (1990) [19]. For our vicious walker model, 
consider the directed graph D = (V, E), where 



V 



= j(ar,i) eZ 2 :x + t = even, t = 0, 1, 2, • • • , t|, 
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and all edges connecting the nearest-neighbor pairs of vertices in V are oriented to the positive direction of t 
axis. Set the weight function 

1 for e= {{x,t- 1) -> (a;- l,t)) 
w(e)={ 1 for e= (\x,t - 1) -» (x + l,t)) 
otherwise. 

Set T e N and Mj = (zj, 0), = (y%,T) e V, i = 1, 2, • • • , TV. Then the Green function of single lattice paths 
from Ui to Vi is 

T 



r i N 

Theorem 5 then gives the Green function of the TV-tuples of nonintersecting lattice paths from I = <^(xt, 0) j 



r 1 N 

to J =\(yi,T)\ as det 



j<N 



(T + Xi - yj )/2 



For x e Z^, and y such that y e Z%, if T e 2N, y + 1 e Z£, if T + 1 e 2N, define 

Mat(T, y|x) = total number of distinct realizations of vicious walk of TV walkers 

from the positions x = (x\, ■ ■ ■ , xn) to y = (j/i, ■ ■ ■ , i/jv) during time T. 

Proposition 4 is now generalized as follows. 



Proposition 6 



AMT, y|x)= det 



T 



(T + Xi - yj )/2 



3 Diffusion Scaling Limit 



Recall that ({S(t)} t=0i i i 2,...,T, Qt) denotes the vicious walk with the noncolliding condition up to time T > 
starting from the positions x e Z< . For L > 1, we consider probability measures /j,* t on the space of continuous 
paths C([0,T] -» K w ) defined by 

where S(i),i > 0, is now considered to be the interpolation of the TV-dimensional random walk S(t),t = 
0,1,2,.... We study the limit of the probability measure ^ T , L — > oo. 
We put 



X\ < X2 < 



< x N 



which can be called the Weyl chamber of type Aw-i (see, for example, [6]). By virtue of the Karlin-McGregor 
formula [10, 11], the transition density function /jv(i,y|x) of the absorbing Brownian motion in R< and 



the probability A/jv(t, x) that the Brownian motion starting from x e M< does not hit the boundary of k< up 
to time t > are given by 



M*,y|x) = det 



i 



/27rf 



-(x 3 -y,) 2 /2t 



, x,yeR^, 



and A/jv(£,x) = / dy/jy(i, y|x), respectively. We put /ijv(x) = jT (x,- — a;,), and let = (0,0, ■■-,0) 



l<i<j<N 

l N , which describes the state that all TV particles are at the origin 0. 
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Theorem 7 (i) For any fixed x £ Z< and T > 0, as L — > oo 7 /i* T (-) converges weakly to the law of 
the temporally inhomogeneous diffusion process X(t) = (Xi(i), X 2 (t), . . . , X N (t)), t € [0, T], with 
transition probability density gN,r{s, x; t, y); 

g N .r(0, 0; t, y) = c^T^ 1 )/ 4 ^ 2 / 2 exp |-^| h N (y)Af N (T t, y), 

, , x f N (t-s,y\x)Af N (T-t,y) 

g N , T (s,x;t,y) = —- , 

Mn(i - s,xj 

for < s < t < T, x, y G R< , w/iere Cjv = 2" JV / 2 / ]T£Li r(i/2) wit/i i/ie gamma function T. 
(ii) T/ie diffusion process X(t) solves the following equati" 



ton: 

N, 



dXi(t) = dBi(t) + bf(t,X(t))dt, te[0,T], i = 1, 2, . . 
w/iere Bi(t), i — 1,2, ... ,N, are independent one- dimensional Brownian motions and 



bf(t, x) = A In A^v(T - t, x), z = 1, 2, . . . , N. 



Figure 4 illustrates the process X(t), i G [0, T], when all N particles start from the origin; X(0) = 0. 
t 




Figure 4: The process X(t), t e [0, T], starting from 0. 

Although Ml.t(') i s tne probability measure defined on C([0,T] — » M w ), it can be regarded as that on 
C([0,oo) -> R w ) concentrated on the set {w e C([0,oo) -> R N ) : tu(i) = w(T), t > T}. Next we consider the 
case that T = T(L) goes to infinity as L — > oo. 



Corollary 8 (i) Lei T(L) be an increasing function of L with T(L) — > oo as L — > oo. For any /ixed x S Z< ; 
as i — > oo 7 /x* T ( L )(-) converges weakly to the law of the temporally homogeneous diffusion process 
Y(i) = (Yi(i), Y 2 (t), . . . , YN(t)),t € [0,oo), wit/i transition probability density pat(s, x; t, y); 



PA r(0, 0; t, y) = c^"^ 2 exp {-^} My) 2 , 
p N {s,x;t,y) = 1 /jv(t - s,y|x)/i JV (y), (3.1) 

«Af(Xj 

/or < s < i < oo, x,y e w/iere = (2 7 r)- Ar / 2 / fljli r(i). 

(ii) T/ie diffusion process Y(t) solves the equations of Dyson's Brownian motion model with the param- 
eter (3 = 2, 

dYi{t) = dBi(t) + Y(t)-Y(t) dt ' * e [°'^)' * = 1,2,...,JV. 
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Here we only give the proof of a key lemma used to prove Theorem 7, in order to demonstrate that the 
Schur function plays an important role. See Katori and Tancmura [12] for the complete proofs of Theorem 7 
and Corollary 8. 

For L > we introduce the following functions: 



H(x) = 2 



Lx 



x e M, and <p L {x) = ($l{x\), <Pl{x2), 4>l{xn)^, 



x e 



where [a] denotes the largest integer not greater than a. Let Vjv(T, y|x) = 2 nt Mn(T, y|x), where Mn(T, y|x) 
is given by Proposition 6. 



Lemma 9 For t > 0, x e Z^ and y e 



V N (<l>L*{t)AL{y) x) 

= c^t-^/ 2 ^ g) cxp {-^} ^(y) (l + 



O 



as L 



Proof. It will be enough to consider the case that x = 2u = (2u\, ■ ■ ■ , 2um) & Z< , y = 2v = (2v\, • • • , 2vn) € Z< , 
and </> L 2 (i) = 2/,/eZ + , where Z+ = {1, 2, 3, • • •}. Then 



and 



M N L L 2(t),<j> L (y) x) = M w (2£,2v|2u) = det 

(2l)\ 



21 

£ + Uj - Vi 



21 

l + Uj — Vi 



{l + Uj - Vi)\{l - Uj + Vi)\ 

(2£)\ 



with 



Aij(£,v,u) = 



{£ + Vi - uj + l) Uj 



Aij(£,v, u), 



{£-Vi + l) Uj ' 

where we have used the Pochhammer symbol; (a)o = 1, (a)i = a(a + 1) • • • (a + i — 1), i > 1. Then 



Aij(e,v,u) 



(3.2) 



The leading term of det 

l<i,j<N 



Aij(e,v,u) 



D\ (v, u) = det 

l<i,j<N 



in L — > oo is 



( _ 1} iV(iV-l)/2 

l<i,j<N 



i + Vj 
t-Vi 



Let £(u) = (£i(u), . . . ,£jv(u)) be a partition specified by the starting point 2u defined by 

&(u) = ujv-i+i — (N — j), j = 1,2, . . . ,N. 

We have 



A(v,u) =(-i)^-D/> det 

1<J,J<JV 



N—j 



a«(u) 



£ — Ul ' ' £ — V]y 
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(-1) 



= n 



N(N-l)/2 



n 



l<i<j<N 

2£(vj - vi) 



+ Vi £ + Vi 



£ — Vi £ — Vi 



S£(u) 



e + vi e + v N 

£ — vi ' ' £ — vn 



(£-Vi)(£-Vj) « (u) 



Vl 



£ - v N J ' 



l<i<j<N 

where we have used Lemma 1 for the Schur function associated to the partition £(u). Proposition 3 gives 

&(u)-£j(u)+j-i 



« e(u) (i,i,...,i)= n 

l<i<j<N 

Therefore the leading term of £>i(v, u) in L — > oo is 

2(vj - Ui) 



D 2 (v,u) 



n " 

l<i<j<N 



x s £(u) (l, !,...,!) 



= r , " f '- 1 »/¥(^ 1 »/ 2 i 1I( (v)i !l »(u) [] — 



l<i<j<N 



J -i 



i—l v y 



On the other hand, by Stirling's formula we see that 
(2£)\ 

i=l 



nTT- 



2\ -<-V2 



1 ~ 
1 + 



L-L (^- Wj )!(^ + v0! 
From (3.2), (3.3) and (3.4) 

Vn (fa (t) , 4>l (y) | x) = 2- 2W M« (fo 2 (t), (y) |x) 



I - O I - 



(3.3) 



(3.4) 



AT/2 



Q)M2u) 



exp < - 



1 + 



= c 



JV 



JV 



2t 1 + °U 



Then we obtain Lemma 9. 



x e K< , then 



Corollary 8 is obtained from Theorem 7 by the following evaluation of asymptotic [12]. Let t > and 

JV 
< ' 

W , X ) . ± hN (^) (, + a (M)) , i„ the limit M _ o, 



where cjv = ^ w/2 IL=i ( r W/ r (V 2 )}- 

4 Eigenvalue Process of Hermitian Matrix- valued Processes 

We consider complex-valued processes £ij(t) € C, 1 < i,j < N,t G [0,oo), with the condition £ji(t)* = 
£ij(t), and introduce Hermitian matrix-valued processes E(t) — [£,ij(t)) . We denote by U(t) = 

V / l<i,j<N 

\Uij(t) ) the family of unitary matrices which diagonalize so that 

V / l<i,j<N 

U(ty~(t)U(t) = A(t) - diag{Ai(t), A 2 (t), • • • , AjvW}, (4.1) 

where |Ai(t)| are eigenvalues of E(t) and we assume their increasing order 

Ai(i) < A 2 (t) < ••• < Ajv(t). 
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Define T t3 (t), l<i,j<N, by 



T ij {t)dt = (u(t)Uz(t)U(t)^^(u{t)Uz(t)U(t)^^ (4.2) 
where dZ(t) = [d^ij(t)\ . The indicator function l^y gives Ij^x = 1 if the condition lu is satished, and 

^ / l<i,j<N 1 



= otherwise. 



Theorem 10 Assume that £ij(t), 1 < i < j < Af, are continuous semimartingales. The process of eigenval- 
ues X(t) = (Ai(i), \2(t), • • • , Ajv(O) satisfies the stochastic differential equations 



d\i(t) = dMi(t) + dJi(t), t e [0,oo), i = l,2,---,N, 
where Mi{t) is the martingale with quadratic variation 



(Mi) t = / T H (s)ds 
Jo 

and Ji (t) is the process with finite variation given by 

N 1 

dJi{t) = ^2 x~(tj - x 3 {t) 1{xdt) ^ {t)}Tl ^ tS)dt + dTt ^ 
where dTi(t) is the finite-variation part of (iJ (t)^ dE(t)U (t)j . 



(4.3) 



(4.4) 



(4.5) 



This theorem is obtained by simple generalization of Theorem 1 in Bru [1]. A key point to derive the theorem is 
applying the Ito rule for differentiating the product of matrix- valued semimartingales: If X and Y are N x N 
matrices with semimartingale elements, then 

d(X^Y) = (dX)^Y + X\dY) + {dX)\dY). 

We give the proof in Appendix B. (See Remark 3 below.) 

Let Bij(t), Bij(t), 1 < i,j < N, be independent one-dimensional Brownian motions. For 1 < i, j < N we 

set 



Sij {t) = < 



■^=B i0 {t), if i<j, 



Bu(t) 
1 



if i = j, and if) = < 



V2 



B^t), if i>j, 



1 

72 



—B i:j (t), if i<j, 



0, 



if i = j, 



I ~V5 Bji ^' [ii>j ' 



A Hermitian matrix-valued process is defined by 

s(t) = Uij(tj)„. =(*«(*) + ^i aij (t)) te[o,oo). 

V / 1<i,j<N \ / 1<i,j<N 



(4.6) 



By definition d£ij(t)d£ki(t) — Sudjkdt, 1 < k,l < N, and thus Tij(t) = 1. Theorem 10 thus implies that the 
eigenvalue process X(t) of this matrix- valued process (4.6) solves the equations of Dyson's Brownian motion 
model with the parameter = 2 [3] 



dXi(t) = dBi(t) 



E 



dt, te [0,oo), i = 1,2, - ■ ■ , AT, 



(4.7) 



where Bi(t), i — 1, 2, • • • , N are independent one-dimensional Brownian motions. 
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Remark 3. In general, Eqs.(4.3) with (4.4) and (4.5) for the eigenvalue process X(t) = (Ai(t), • • • , A;v(t)) 
depend on unitary matrix U(t) through Tij(t) defined by (4.2). The equations written in the form, 

d\i(t) = J2 a ij(^ x (t))dBj(t) + Pi(t, X(t))dt, t e [0, oo), i = 1, 2, • • • , N, (4.8) 

3 

where the coefficients aij(t,X) and (3i(t,X) are functions not only of A but also of other variables, are 
generally called stochastic differential equations (SDE's) in [9] (see Definition 1.1 with Eqs.(l.l), 
(1.1') in Chapter IV 'Stochastic Differential Equations' on page 159.) In the special case, in which these 
coefficients are only depending on X(t), equations are given in the form 

dXi(t) = J2 a v(Mt))dB 3 (t) +bi(X(t))dt, t€ [0,oo), i = 1,2, •••,7V, (4.9) 

3 

and they are said to be of the Markovian type (see page 172 with Eq. (2.11) in [9]). The condition that 
the SDE's of eigenvalue process are reduced to be of the Markovian type may be that the matrix- valued 
process S(i) is unitary invariant in distribution. By virtue of properties of Brownian motions, the 
Hermitian matrix-valued process S(t) defined by (4.6) is unitary invariant in distribution, and thus the 
obtained SDE's of Dyson's Brownian motion model are of the Markovian type. 



5 Concluding Remarks 

Corollary 8(ii) and Eq.(4.7) with (3 = 2 implies that the temporally homogeneous process Y(t) obtained as a 
diffusion scaling limit of vicious walks and the eigenvalue process X(t) of the Hermitian matrix-valued process 
(4.6) are equivalent in distribution. The formula (3.1) in Corollary 8 shows that it is the /i-transform in 
the sense of Doob [2] of the absorbing Brownian motion in the Weyl chamber R<, since /ijv(x) is a 
strictly positive harmonic function in R< [7] . An interesting relationship between this temporally homogeneous 
process Y(t) (Dyson's Brownian motion model with the parameter (3 = 2) and the temporally inhomogeneous 
process X(t) given by Theorem 7 was reported in [12, 14]. A systematic study on the relations among various 
matrix- valued processes, standard, chiral and non-standard random matrix theories, and families of noncolliding 
diffusion processes was reported in [15]. 

For the noncolliding diffusion processes starting from 0, the multi-time correlation functions were 
calculated by Nagao and the present authors using the quaternion determinants of self-dual quaternion matrices 
(i.e. pfaffians) and the scaling limits of the infinite particles N — > oo and the infinite time- interval T —> oo 
were investigated [18, 13]. Further study of infinite systems of noncolliding diffusion particles will be 
reported elsewhere [16]. 



A Proof of Theorem 5 

By definition of determinant 



dct 

l<i,j<N 



G(ui,Vj) = ^2 s g n ( CT ) G ( u i;Mi))''' G ( M ^'^w)> ( A -!) 

<tG6jv 



where ©at is the set of all permutations of {1, 2, • • • , TV}. We may interpret (A.l) as a generating function for 
(N + l)-tuples {a,P u ---,P N ), where a € & N ,Pi & V{u u <v (i) ), i = 1,2,---,A. 

Consider an arbitrary configuration (a, P\, ■ ■ ■ , P/v) with at least one pair of intersecting lattice paths. We 
set the order of all vertices in V. Let v denote the last vertex among all vertices that occur as points of 
intersection among the lattice paths. Among the lattice paths that pass through v, assume that Pi and Pj are 
the two whose indices i and j are the smallest (see Figure 5). 

Write 

P = P (-» v)Pi(v ->), and Pj = Pj (-» v)Pj («-»). 
For a configuration (a, Pi, ■ ■ ■ , Pn), define as shown in Figure 6 

Pi = p(^ v )Pj( v ^), 
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Figure 6: (a) Paths P t and P r (b) Paths P[ and Pj. 

Figure 5: A configuration (a, P\, ■ ■ ■ , P/v). The last in- 
tersecting vertex is denoted by v and lattice paths Pi 
and Pj are chosen, both of which pass through v. 

p; = p (^ v)p t { v 

P' k = P k fork ^1,3, 

a 1 = cro (i, j), where denotes an exchange of i and j. 

The operation 

(a,P 1 ,---,P N )^(a',P{,---,P' N ) 

preserves the set of vertices of intersection, and is an involution. The weight of lattice paths is the same, but 
the sign is changed. So any such pair j(er, Pi, • • • , Pjv), {&', P{, • • • , Pjv)} a PP ear in (A.l) is canceled out. 

The only configurations remain in (A.l) arc nonintcrsccting lattice paths. Since I is assumed to be in- 
compatible with J, for nonintersecting lattice paths a — id, i.e., sgn(cr) = sgn(id) = 1. I 

B Proof of Theorem 10 

We consider a matrix-valued process A(t) = (a^ (t))i<i.j<Ar defined by 

dA{t) = U(t)Uu(t) + ^dU(t)UU(t), t e [0, oo) 

with A(0) = 0. Since U(t)^U(t) — In for all i, where In denotes the N x N unit matrix, 

= d(U(t^U(t)) = dU(tfU(t) + U{t)UU(t) + dU(t)Uu(t). 

Then 

dA(t)t = dU(tfU{t) + ^dU(t)Uu(t) 

= -U{t)UU(t) - ^dU(t)UU(t) = -dA(t), 

that is, dA(t) is anti-Hermitian. We also see that 
-dA(t)dA(t) = dA{t)UA(t) 

u(t)Uu(t) + l-du(t)Uu(t)] (u(t)Uu(t) + \dU{t)Uu{t) 



2 

= dU{tfU(t)U{t)UU{t) = dU{t)UU{t). (B.l) 
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This implies 

dU(t) = U(t) (dA(t) + ^dA(t)dA(t)J . (B.2) 

By (4.1) 

dA(t) = dU{tfZ(t)U(t) + U(t)Uz(t)U(t) + U(tfZ(t)dU(t) 

+dU(t)UZ(t)U(t) + dU(tfZ(t)dU(t) + U(t)Uz(t)dU(t) 

= U(tfdZ{t)U{t) + ^A(t)U(t)Uu(t) + (A(t)U(t)Uu(t))^ 

+I^U(t)Uz(t)dU(t) + {U(t)Uz(t)dU(t))^ +dU(t)' i E{t)dU(t). 

Each term in the RHS is rewritten as follows: 

A(t)U(tfdU{t) = A(i) (u(tfd(7(t) + ^dU(t)UU(tyj - ^A(t)dU(t)Uu(t) 



= A{t)dA{t) + ^A(t)dA(t)dA(t), 



where (B.l) was used, and 



U(t)Uz(t)dU(t) = U(t)Uz(t)U(t)dA(t), 

dU(tfE(t)dU(t) = dA(t) ] U{tf T ~{t)U{t)dA(t)=dA(tfA(t)dA{t), 



where (B.2) was used. Then we have the equality 

dA(t) = U(tfdZ(t)U(t) + A(t)dA(t) + (A(t)dA(i))t 
+^A(t)dA(t)dA(t) + ^(A(t)dA(t)dA(t)Y 

+U(t)UZ(t)U(t)dA(t) + (U(t)Uz(t)U(t)dA(t)^ + dA{t)^A{t)dA{t). (B.3) 
The diagonal elements of (B.3) give 

d\i(t) = y^uki{t)*uei{t) d£,kt (t) 
k,i 

+2A i (t)d7«(t) + #«(t)+d0«(t)* +dipii{t), l<i<N, (B.4) 
and the off-diagonal elements of (B.3) give 

= u kz{t)*u e: j(t)d£ kl ,(t) + \i{t)daij{t) + \j(t)da 0l {t)* 

k,i 

+Aj(t)d7« (*) + x j (*)*>&(*) + d<j>a (t) + + dipij (t), 1 < i < j < N, (B.5) 

where we have used the notations 

diijit) = (^dA(t)dA(t^j =^J2da ik (t)da kj (t) = d^t)*, 

#«(t) = (u(tfcE{t)U(t)dA(tj) = Uki(tyd^ M (t)u em (t)da mj (t), (B.6) 

^ k,£,m 

#«(t) = (dA(tf A(t)dA(t)y 

= ^da ki (t)*\ k (t)da kj (t) = -^da tk (t)\ k (t)da kj (t). 
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Since {'Jij(t))i<i,j<N , {4'ij{t))i<i,j<N , (i>ij{t))i<i,j<N are functions of finite variations, (B.3) gives 

= (u(t)Uz(t)u(t)) (u(t)Uz(t)u(t)) 



d{Mi) 



= ^2^2 u kt {i)*u u (t)u mi (t)*u ni (t)d£ ke {t)d$ mn (t) . 

k,l m,n 



On the other hand, for dA(t) is anti-Hermitian, (B.5) gives 

(Aj (t) - A, (t))da y (t) = 51 u « (*) 



+ (Ai(t) + \j(t))d-yij(t) + d<f>ij{t) + d(j)ji{t)* + #«(t), (B.7) 



and using this equality we can rewrite (B.6) as 

= y^(Afcft) - \i(t))da ik (t)da k j{t). 
k 

Then the finite-variation part of (B.4) is written as 

2\i(t)dj u (t) + d(j) U (t) + d(pu(t)* + dipu(t) 

= E j A *W + 2 ( A ^*) - A »W) - (*) | (t)daji (t) 
= J2^j(t) - XiiWaij^dajiit) 



= Yl \.( t ) 1 x-(t) 1{XiW¥:Xj(t)} ^ u k t {t)*u e:j (t)u m: j(tyu m (t)d£_ k i{t)d£_ mn {t) 7 

j ^ ' ' k,e,m,n 

where (B.7) was used in the last equation. This completes the proof. 
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